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ABSTRACT 



A full-sky template map of the Galactic free-free foreground emission component 
is increasingly important for high sensitivity CMB experiments. We use the recently 
published Ha data of both the northern and southern skies as the basis for such a 
template. 

The first step is to correct the Ha maps for dust absorption using the 100 /im dust 
maps of Schlegel, Finkbeiner & Davis (1998). We show that for a range of longitudes, 
the Galactic latitude distribution of absorption suggests that it is 33 per cent of the 
full extragalactic absorption. A reliable absorption-corrected Ha map can be produced 
for ~ 95 per cent of the sky; the area for which a template cannot be recovered is the 
Galactic plane area |6| < 5°, Z = 260° — 0° — 160° and some isolated dense dust clouds 
at intermediate latitudes. 

The second step is to convert the dust-corrected Ha data into a predicted radio 
surface brightness. The free- free emission formula is revised to give an accurate ex- 
pression (1 per cent) for the radio emission covering the frequency range 100 MHz to 
100 GHz and the electron temperature range 3000 to 20000 K. The main uncertainty 
when applying this expression is the variation of electron temperature across the sky. 
The emission formula is verified in several extended Hii regions using data in the range 
408 to 2326 MHz. 

A full-sky free-free template map is presented at 30 GHz; the scaling to other 
frequencies is given. The Haslam et al. all-sky 408 MHz map of the sky can be cor- 
rected for this free- free component, which amounts to a « 6 per cent correction at 
intermediate and high latitudes, to provide a pure synchrotron all-sky template. The 
implications for CMB experiments are discussed. 

Key words: cosmic microwave background - radio continmmr: ISM - dust, extinction 
- Hii regions - radiation mechanisms: thermal 



1 INTRODUCTION 

Current CMB experiments are sensitive enough to measure 
the primordial fluctuations which have an amplitude in the 
range of ~ 20—100 over the £-range 10 — 2000 which cor- 
responds to angular scales of ~ 10° — 10 arcmin (Hanany et 
al. 2000; Mauskopf et al. 2000; Padin et al. 2001; Halverson 
et al. 2002; Scott et al. 2002). The angular power spectrum 
of these fluctuations contains a wealth of cosmological infor- 
mation. One of the crucial factors in determining an accurate 
power spectrum of these fluctuations is understanding and 
removing the foreground contamination. The amplitude of 
the foreground signal depends on frequency, angular scale 
and region of sky being observed. CMB foregrounds com- 
prise point sources and diffuse Galactic foregrounds. Point 
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sources are a particular problem at smaller angular scales 
(< 30 arcmin) , while diffuse foreground structures become 
dominant on larger angular scales. Here we concentrate on 
the problem of diffuse Galactic foregrounds at frequencies 
below < 100 GHz, a frequency range which is used in cur- 
rent and upcoming CMB experiments such as the MAP and 
Planck satellites. 

1.1 Galactic foregrounds 

The diffuse Galactic foreground has 3 (possibly 4) com- 
ponents, i) Synchrotron emission from relativistic electrons 
spiralling in the Galactic magnetic field dominates at fre- 
quencies below 1 GHz with a spectral index (T oc u^^) 
of /3 « 2.7 - 3.2 (Davies et al. 1996). It is traced by low- 
frequency surveys such as those of Haslam et al. (1982), 
Reich & Reich. (1988) and Jonas et al. (1998) where it is 
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seen to extend well above the plane - e.g. the North Polar 
Spur which extends to Galactic latitude b « 80° . ii) Free-free 
bremsstrahlung emission from thermal electrons has a flatter 
spectral index of /3 ~ 2.1. It is the dominant foreground at 
frequencies between u — 10 — 100 GHz, where microwave 
CMB experiments operate (e.g. DASI, CBI, VSA, MAP, 
Planck). The optical Ha line is a good tracer of free- free 
emission although it requires corrections for dust absorption, 
iii) The vibrational emission from thermal dust is dominant 
at frequencies above ~ 100 GHz. This is thermal emission 
from warm dust and is well-traced at A ~ 100 /xm where 
it has its peak, iv) An anomalous component has recently 
been discovered (Kogut et al. 1996; Leitch et al. 1997; de 
Oliveira-Costa et al. 1997,1998,1999,2000,2002; Finkbeiner 
et al. 2002). Draine & Lazarian (1998) proposed that this 
could be due to spinning dust grains emitting in a 1-2 oc- 
tave band centred at ~ 20 GHz. It was flrst believed to be 
due to free- free emission (Kogut 1996), but this has been 
ruled out primarily by the lack of associated Ha emission. 

The separation and quantifying of the individual CMB 
foreground components is a continuing challenge. Each com- 
ponent is of interest in its own right; its angular distribution 
and spectrum are basic parameters. They all contribute as 
foregrounds to the CMB and their removal is necessary to 
realize the full potential of current high sensitivity CMB sur- 
veys where Ci values are to be measured to an accuracy of a 
few percent. Such separations are best made using templates 
for each component. With the advent of the Ha survey, such 
a template is now available for each of the 4 proposed com- 
ponents. Although some approaches to separation of compo- 
nents are being made which reduce the assumptions about 
the template parameters (e.g. Baccigalupi et al. 2000; Vielva 
et al. 2001) they still face difficulties where the components 
are correlated. For example all the foregrounds are quasi- 
correlated as they increase strongly towards the Galactic 
plane. Also, the Ha and (spinning) dust emissions show con- 
siderable correlation over most of the sky (see Section |7|l; 
this was at the root of the original misunderstanding of the 
anomalous component. If one has a good free-free template 
as we derive here, the other 3 foreground components can 
be well-characterized except possibly at the lowest Galactic 
latitudes. 

1.2 Large-area Ha surveys 

Extensive filamentary Ha emission regions, of galactic ori- 
gin, were first found by Meaburn (1965, 1967) extending 
to high Galactic latitudes; Sivan (1974) demonstrated the 
wealth of Ha emission at low-intermediate latitudes but at 
a relatively low sensitivity of « 15 R; we note that 1 Rayleigh 
{R) = 10'^/47r photons s'^ cm'^ sr'^ = 2.41 x 10^ erg s'^ 
cm"^ sr"^ = 2.25 cm"®pc for Te = 8000 K gas. Several large 
and very sensitive Ha surveys are now well under way. The 
most relevant surveys are listed in Table Q In the northern 
sky covering declinations 5 ^ —30° the Wisconsin H Alpha 
Mapper project (WHAM) is the most sensitive Ha survey 
to-date, but is restricted to an angular resolution of ~ 1° 
(Reynolds et al. 1998; Haffner 1999). The final sensitivity 
in each field is « 0.05 R and the data are calibrated to an 
accuracy conservatively estimated at 10 per cent. WHAM is 
a dual etalon Fabry-Perot spectrometer with a velocity res- 
olution of 12 km s"'^. The spectra can be used to remove the 



geocoronal Ha emission from the Earth's upper atmosphere 
which varies in strength from 2-13 R; it is much brighter 
than the Galactic signals at high Galactic latitudes. 

The ongoing Virginia Tech Spectral-line Survey (VTSS) 
is a complimentary narrow-band imaging survey of the 
northern sky covering 5 ^ —15° (Dennison, Simonetti & 
Topsana 1998). It has a resolution of 1.6 arcmin in each 
10° diameter field, which will be good enough for almost all 
CMB experiments as most of the cosmological information 
is contained on scales between 10° — 10 arcmin. However, 
on larger angular scales, it will most likely be limited by 
the geocoronal emission which appears as a time-varying 
background signal of unknown level. Star residuals which 
have not been subtracted correctly can also be a problem. 
In combination with the WHAM survey, this survey will be 
a powerful tool. 

The recently published (Gaustad et al. 2001) Southern 
H Alpha Sky Survey Atlas (SHASSA) covers the southern 
sky 8 ^ -1-15° at an angular resolution of 0.8 arcmin in each 
13° field. The sensitivity reaches 2 R limited by geocoro- 
nal emission and star residuals which have been partially 
removed by median filtering. Smoothing this survey to a 
resolution of 4 arcmin allows features of about 0.5 ii to be 
detected. The intensity calibration is accurate to about 9 
per cent. 

Two other Ha surveys are in progress. The Manch- 
ester Wide-Field Camera (MWFC) can observe Ha with 
a 32° field of view with 7 arcmin resolution (e.g. Boumis 
et al. 2001) and a sensitivity oi k. 1 R. This instrument is 
particularly sensitive to large scale features (> 1°) which 
may be missed by cameras with smaller fields. We have 
recently observed several selected fields to look for diffuse 
Ha emission, both in the northern and southern hemi- 
spheres. The AAO/UK Schmidt Ha survey has high angular 
resolution (1 arcsec) but is restricted to the Galactic plane 
at relatively low sensitivity (Parker & Phillips 1998) as given 
in Table □ 



1.3 This paper 

In this paper we describe the steps taken to derive a free- 
free emission template from recently published Ha surveys. 
In Section |21 we determine the absorption of the Galac- 
tic Ha emission by dust based on the 100 /xm dust tem- 
plate given by Schlegel, Finkbeiner & Davis (1998), hereafter 
SFD98; a statistical estimate is made of the relative distri- 
bution of the Ha emission and the dust in the line of sight. 
Section |3] summarizes the expressions required to convert 
the dust-corrected Ha emission into microwave emission at 
frequencies relevant to current CMB experiments. The un- 
certainties in the relation are evaluated. This relationship is 
tested in Section |1] using available radio surveys. The final 
free-free template is presented in Section |3 The 408 MHz 
all-sky map of Haslam et al. (1982) corrected for free-free 
emission is derived in Section |S| it is the best available syn- 
chrotron template. We discuss the implications of the new 
template for upcoming CMB experiments and for interstel- 
lar medium (ISM) studies in Section|7| The final conclusions 
of this work are given in Section |S| 
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Table 1. Current Ha surveys relevant to CMB observations 



Survey 


Type 


Field of View 
(degrees) 


Resolution 
(arcmin) 


Sensitivity/ pixel 
(R) 


Sky 
coverage 


WHAM 


Fabry- Perot 
Spectroscopy 


1 


60 


~ 0.05 


S ^ —30° 


VTSS 


Narrow-band 
Imaging 


10 


1.6 


1 


5 ^ -15° 


SHASSA 


Narrow-band 
Imaging 


13.6 


0.8 


KS 2 


S ^ +15° 


MWFC 


Narrow-band 
Imaging 


32 


7 


^ 1 


Selected 
fields 


AAO/ 
Schmidt 


Narrow-band 
Imaging 


5.5 


0.02 


Si 5 - 10 


\b\ ^ 10° 



2 CORRECTION OF Ha FOR DUST 
ABSORPTION 

2.1 The Galactic distribution of absorption 

It is clear that the use of Ha surveys as a template for free- 
free emission will be jeopardized in the Galactic plane where 
visual absorption is typically 1 magnitude per kiloparsec in 
the local arms and where the total absorption to the Galac- 
tic centre is « 20 magnitudes. At intermediate latitudes the 
absorption broadly follows a cosecant law with a vertical slab 
half-thickness of 0.1 to 0.2 blue magnitudes. Such a cosecant 
law was widely used in extragalactic astronomy. In order to 
take account of the known structure in the obscuration, the 
line integral of Hi combined with a factor derived from the 
Shane- Wirtanen galaxy counts was introduced by Burstein 
& Heiles (1978). This procedure took account of the possi- 
bility that the gas-to-dust ratio might not be constant and 
the fact that not all hydrogen is in the neutral atomic form. 
The Burstein-Heiles approach is limited by the angular res- 
olution of Hi all-sky surveys which is currently ~ 30 arcmin. 

A new approach which provides a resolution of 6.1 
arcmin and has a complete all-sky coverage is offered by 
SFD98 who use far infra-red data from the COBE-DIRBE 
and IRAS satellites. This combination of data at a range of 
FIR wavelengths enabled a good zero level to be established, 
an adequate removal of zodiacal hght and a more effective 
destriping than previously. Also, by comparing the DIRBE 
240 /im and IRAS 100 fim to derive a dust temperature, 
an estimate could be made of the dust column density, , 
measured in units of MJy sr~^, in terms of a 100 pm surface 
brightness at a fixed temperature of 18.2 K. This correction 
to a fixed temperature can only be made on the angular 
scale of the DIRBE observations, namely 0°.7, and conse- 
quently is not strictly true on the 6.1 arcmin scale of the 

dust maps. The dust temperatures are typically in the 
range 17 — 21 K which corresponds to a correction of up to 
a factor of ~ 5 in the dust column density (see SFD98). 

2.2 Derivation of a Ha absorption template 

We will adopt the dust template as the indicator of 
dust absorption. This then requires a conversion factor to 



estimate the absorption at 656.3 nm, the Ha wavelength. 
SFD98 use the (B — V) colours of some 470 galaxies at a wide 
range of Galactic coordinates to derive a best fit correlation 
with the value at the position of each galaxy. They find 
the E{B — V) colour can be estimated from with the 
expression E{B - V) = (0.0184 ± 0.0014)D^ magnitudes 
and claim that the reddening estimated in this way has a 
standard deviation of 16 per cent at any given position. 

In order to estimate the absorption at the 
Ha wavelength we use the parametric extinction curve for 
optical wavelengths given by O'Donnell (1994) and find 
the Galactic absorption at Ha to be y4(Ha) — 0.81 A{V). 
Assuming the dust is characterized by a reddening value 
Rv = A{V)/ E{B ~V) — 3.1, this leads to an absorption at 
Ha of 

yl(Ha) = 2.51E{B - V) ^ (0.0462 ± 0.0035)-D'^mag (1) 

where is in units MJy sr~^. The values given in table 6 
of SFD98 are 6 per cent higher than this, corresponding to 
a higher value of Rv- 

A range of reddening laws, characterized by diflerent 
values of Rv, are found in different directions in the Galaxy. 
The different laws are thought to derive from different grain 
chemistry and size distributions. Values of Rv range from 
2.5 to 5, with the higher values often in denser dust clouds. A 
value of 3.1 is widely adopted as representative of the diffuse 
ISM in directions away from dense dust clouds. However, it 
should be emphasized that this is only an average value and 
that any line of sight will have a scatter about this value. 
Putting all these statistical factors together we consider that 
the relationship between yl(Ha) and is accurate to about 
20 per cent. 

2.3 The absorption of Galactic Ha 

The expression given in the preceding section applies to 
absorption of extragalactic objects at optical wavelengths. 
Here we are considering Galactic Ha emission which is mixed 
with the dust. At first sight it might be argued that since 
the dust and gas are uniformly mixed, then the ionized gas 
and the dust would also be uniformly mixed. The correlation 
between gas and dust will be discussed further in Section |7| 
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The major part of the Hq seen from our position in the 
Galaxy is the resuh of ionization by the UV radiation field 
from young stars embedded in the dust plus the contribu- 
tion from nearby regions of localized recent star formation. 
Such regions include the Gould Belt system which reaches 
to Galactic latitudes of 40° or more. 

We seek a first-order absorption formula for correction 
of the Hof emission based on the dust distribution given by 
SFD98. If the ionized gas were uniformly mixed with the 
dust then the average absorption would be half the extra- 
galactic value given in Section 12.21 on the assumption that 
the dust absorption is optically thin (say ^ 0.5 mag). We 
define fd, the effective dust fraction in the line of sight actu- 
ally absorbing the Ha so that the actual absorption of Ha is 
given by fd x ^(Ha), where yl(Ha) is the full extragalactic 
absorption. For example, uniformly mixing of ionized gas 
and dust corresponds to fd ~ 0.5. There would be a sys- 
tematic modification of this value if the z-distribution of the 
dust and the ionized gas were not the same. An example 
would be the ionization produced by the interstellar radia- 
tion field propagating away from the plane and ionizing the 
"under" -side of the gas clouds which were themselves opti- 
cally thick to the UV ionizing radiation. This would lead to 
a narrower z-distribution of Ha emission than dust, which 
is assumed to be well-mixed with the gas, and to a reduced 
value of fd- 

We are now in the position to test the uniform mixing 
hypothesis by comparing at low and intermediate latitudes 
the z-distribution of neutral gas (Hi), ionized gas (Ha) and 
dust. 

2.4 Galactic latitude scans 

By plotting the Ha dependence on Galactic latitude b, and 
comparing this with the unabsorbed gas indicators of Hi 
and FIR dust (D^), we can estimate the Ha absorption as 
a function of b and determine a best-fitting value of fd, the 
effective dust fraction producing absorption. 

For neutral Hi gas we use the 1420 MHz Leiden- 
Dwingeloo northern sky survey (Hartmann & Burton 1997) , 
Ha data from the WHAM survey and the SFD98 100 ^m 

map to trace the dust. Each map was resampled onto a 
Galactic coordinate grid and smoothed to a resolution of 1°. 
Along each l°-wide latitude strip, the values are averaged 
over a range of longitudes to give a representative value for 
that latitude. Fig. shows the latitude dependence of the 
gas and dust in the region I = 30° —60°. Both the Hi gas and 
100 ^m dust follow similar cosecant-law trends. Assuming a 
homogeneous slab of material viewed at different latitudes 
6, the column density along a line of sight is expected to 
follow a cosecant law when viewed from the central plane; 

I ^ Ao+Ai/sm{\b\) (2) 

where Ao is an offset and Ai is the amplitude of the cosecant 
law. This is indeed found to be the case, with deviations oc- 
curring due to large concentrations of emitting material such 
as extended Hll regions and features like the Local System 
(Gould Belt) . The region I — 30° —60° is one of the "cleaner" 
regions of the plane where no large structures exist. This re- 
gion includes are own spiral arm and the Sagittarius arm; 
here the dust emission on the Galactic ridge reaches an av- 
erage of 500 MJy sr~^ corresponding to A{lla)= 23 mag for 



an extragalactic object. At intermediate and high latitudes 
(|6| > 15°), is < 8 MJy sr~^ and the Ha absorption is 
therefore small. 

The best-fitting cosecant law is a good match to the Hi 
and dust distributions as shown in Fig. Q This figure also 
shows that the Ha intensity is greatly attenuated at low 
latitudes (|6| < 15°) due to the increased dust absorption 
nearer to the plane. Note also that there is a significant off- 
set ^o; this corresponds to the Sun being in a local "hole". 
The Ao/Ai values suggests that there is a ~ 30 — 50 per cent 
deficit in the local slab density. If one adopts equation Q as 
the model for Ha absorption we can correct the Ha emission, 
adopting a value of fd as discussed in Section [2.31 We can 
then determine a value for fd by correcting the Ha latitude 
distribution for absorption, so that it matches that of the 
gas and dust distributions shown in Fig Q By varying fd 
we find the best-fitting value is given by fd — 0.331q jg. 
The fit was made by deriving the cosecant law using data 
at |ti| > 20° and fitting fd to this cosecant law over |fe| = 5° 
to 15°. For |6| < 5°, the absorption is likely to be > 1 mag 
and therefore any derived value for fd becomes unreliable. 
The value fd = 0.33 derived here is representative of the so- 
lar neighbourhood of our Galaxy within a few kilo-parsecs. 
Other longitudes gave similar results, although the cosecant 
fitting is difficult due to local structures such as the Gould 
Belt system and other bright extended Hll regions. Further 
information on the value of fd is obtained from the study 
of the radio free-free emission from several nearby Hll com- 
plexes given in Section |1] 

The value of fd derived in this section is not consistent 
with uniform mixing of ionized gas and dust which requires 
fd = 0.5. However, there is a large uncertainty in fd and 
it is expected that fd will vary significantly across the sky. 
However, as a first-order approach, we will adopt the value 
fd = 0.33 when estimating the free-free template in Section 
|5] This lower value is indicative of non-uniform ionization; 
for example, a region of cloud is ionized by the interstellar 
radiation field which is in a narrower Galactic layer of 0,B 
stars. 



2.5 The Ha absorption template 

Using equation Q we can convert the dust map to an 
absorption map which can be applied to the Ha intensity 
/hq, using 

rcorrcctcd r w 1 n-D'^ X 0.0185 X /o\ 
lua = J Ha X 10 " (3) 

The value of fd is taken to be 0.33. The absorption template 
to be applied to Galactic Ha and smoothed to a resolution 
of 1° is shown in Fig.|21 

The template clearly shows that the absorption of 
Ha by dust is modest - less than 0.2 mag (20 per cent cor- 
rection) for most of the sky. Beyond \b\ > 50°, it is less than 
0.05 mag (5 per cent correction). At |6| < 5°, the absorp- 
tion increases beyond 1 mag for much of the Galactic plane. 
Individual dust clouds can have > 1 mag of absorption at 
intermediate latitudes. For a major portion of the sky the 
dust correction clearly needs to be taken into account. For 
example, the correction of more than 50 per cent is required 
in a limited area of the sky such as the Gould Belt system. 
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Figure 1 . Galactic latitude scans for the longitude range I = 30° — 60° . (a) Hi from Leiden-Dwingeloo survey data, (b) Dust emissivity at 

100 /im from the SFD98 template, (c) Ho data from the WHAM survey. In (c) the Ho data corrected for dust absorption {f^i = 0.33) 
are shown as filled circles while the uncorrected data are open circles. The best-fitting cosecant law is shown for each distribution (see 
text). Note the logajrithmic scale on the vertical axes. 
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Figure 2. The absorption template to be applied to the Galactic H« emission. This assumes an absorbing dust-fraction of = 0.33, as 
derived in Section 12.41 and the the total extragalactic value A{Ra) given in equation 0. Linear grey-scale and contours show absorption 
in magnitudes. Contours are given at 0.05 (dot-dashed), 0.1, 0.2, 0.3, 0.4, 0.5 and 0.7 mag. Absorption above 1 mag is black on the grey-scale. 



It is interesting to note that the absorption can be low 
(< 0.5 mag) close to the Galactic plane in the longitude 
range 160° < I < 260°. These regions may be useful in test- 
ing the dust absorption and Ha templates at low latitudes. 
However, for the remainder of the Galactic plane < 5°) 
the absorption is too great to make any reasonable estimate 
of Ha emission. 



The uncertainty in deriving the Hq absorption from 
the dust using equation Q is considered to be 20 per cent 
(see Section [2.211 . A second factor in the uncertainty is the 
relative distribution of gas (Hq) and dust in the line of 
sight (/d). Near the Galactic plane where there are many 
clouds/features in the line of sight, our value of fd = 0.33 
is realistic, while at higher latitudes where there are only a 
few clouds in the line of sight, the uncertainty (in a smaller 
total absorption) is greater. If we adopt a 20 per cent un- 
certainty in fd and combine it with a 20 per cent error in 
equation we obtain a total uncertainty of 30 per cent in 
the correction to the Ha intensity for dust absorption. Thus 
if there is an Ha absorption ( fd x ^(Ha)) of 1.0 mag, the 30 
per cent uncertainty in multiplying factor required to give 
the estimate of the true Ha intensity is 2.5 ± 0.7. However 
for a total absorption of 0.2 mag the multiplying factor is 
1.20 ± 0.07. This illustrates the higher fractional accuracy 
of estimating the Ha intensity at lower levels of absorption. 
In summary, the absorption of Ha can be estimated up to 
a limit of ~ 1 mag; this corresponds to = 65 MJy sr~^ 
for fd ~ 0.33 and is shown as black on Fig. |5| 



3 CONVERSION OF Ha TO FREE-FREE 
CONTINUUM EMISSION 

The ionized interstellar medium generates both radio free- 
free continuum and Ha emission. Ionized hydrogen alone 
is responsible for the Ha emission while ionized hydrogen 
and helium produce the radio continuum. Both the optical 
and radio emission are functions of electron temperature, 
Te- A relationship between the two can be determined on 
certain assumptions. We will derive this relationship and 
give an estimate of the uncertainties which come from the 
theory and from the variation in the electron temperature 
at different locations within the Galaxy. 



3.1 Ha emission 

The Balmer line emission from ionized interstellar gas is well 
understood; a clear discussion is given by Osterbrock (1989). 
However, the Ha line intensity depends upon whether the 
emitting medium is optically thick (case B) or optically 
thin (case A) to the ionizing Lyman continuum. In the Hll 
regions and nebulae studied using Ha it is believed that 
case B applies (Osterbrock 1989). Knowing the emission 
coefficients for both case A and case B for H/3 and us- 
ing the Balmer decrements (Hummer & Storey 1987), the 
Ha emission can be readily calculated for the tempera- 
ture range 5000 — 20000 K both for case A and case B. 
Fig. 13 shows the Ha intensity per unit Emission Measure 
(EM = J nidi) for Te = 5000, 10000, 20000 K for both case 
A and case B. 

The emission coefficients are generally accepted to be 
accurate to ~ 1 per cent (Pengelly 1964; Hummer & Storey 
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Figure 3. Ha emission per unit Emission Measure as a func- 
tion of Te for case A (optically thin to Lyman photons) and case 
B (optically thick to Lyman photons) for = 10^ cm~'^. The 
circles represent the values calculated by theory given by Hum- 
mer & Storey (1987). The curves are the expressions given by 
Valls-Gabaud (1998); case B is given in equation Bl . 

1987). For case B, the intensity varies very weakly with num- 
ber density, amounting to a few percent over 2 orders of 
magnitude in density. The values shown in Fig. |H] are for 
a density (rie = lO'^ cm~^) gas. Case A gives « 30 per 
cent lower intensity than case B over the hkely temperature 
range. 

The situation for intermediate and high latitudes ap- 
pears to favour case B on the following arguments. For case 
B to apply, the optical depth r, to Lyman continuum pho- 
tons (A J; 900 A) must be greater than unity in the emit- 
ting region. Pottasch (1983) gives the mean free path of 
Lyman continuum photons as ~ 0.052ng^ pc, correspond- 
ing to an optical depth relation r = 19.2nHl where / is path 
length in pc. For a typical Hll region in the Galactic plane 
He ~ riH = 10^ — 10® cm~'^ on typical scales of 1 — 10 pc and 
case B conditions are clearly satisfied with r ~ 10'' — lO^''. 
At intermediate latitudes faint Ha features with EAI ~ 1 
cm~® pc can be identified with scale sizes of ~ 10 — 100 pc. 
The corresponding Ue is 0.1 — 0.3 cm~^. The optical depth 
to Lyman continuum in these low density regions is then 
1 — 30, still satisfying the case B conditions. 

A useful and accurate expression for the Hq intensity 
based on data from Hummer & Storey (1987) is 
given by Valls-Gabaud (1998) for case B in units of 
erg cm~^ s~^ sr~^: 

/(Ha) °"=i?^ 9.41 X io-8r-i oi^lO-" °23/^^ x (i?M),^-6p,(4) 

where T4 is the electron temperature in units of 10^ K. This 
is plotted as a continuous curve along with the equivalent 
curve for case A in Fig. |3] which shows that the agreement 
of equation @ with theory is better than 1 per cent over 
the temperature range 5000 - 20000 K. 

3.2 Radio continuum emission 

The free-free radio continuum emission for an ionized gas at 
Te < 550, 000 K and in Local Thermal Equilibrium (LTE) 
is well described in terms of a volume emissivity in units of 
erg cm"^s"^Hz"^ (Oster 1961): 



Table 2. Gaunt factors for relevant frequencies and electron tem- 
peratures calculated from Hummer (1988) and accurate to 0.7 per 
cent. 



u (GHz) 


4000 


6000 


(K) 
7000 


8000 


10000 


0.4 


5.82 


6.03 


6.11 


6.18 


6.39 


1.4 


5.16 


5.37 


5.45 


5.52 


5.73 


2.3 


4.90 


5.11 


5.19 


5.26 


5.47 


10 


4.13 


4.34 


4.42 


4.49 


4.70 


30 


3.55 


3.76 


3.84 


3.91 


4.12 


44 


3.35 


3.56 


3.64 


3.71 


3.92 


70 


3.10 


3.32 


3.40 


3.45 


3.67 


100 


2.92 


3.13 


3.21 


3.28 


3.49 



e^-' = 6.82 X 10 Z ncUionT^ e ' < gff > (5) 

Here < gff{Z,v,Ts) > is the velocity-averaged Gaunt fac- 
tor. The Gaunt factor comprises all the terms by which 
the quantum mechanical expressions differ from the classical 
ones. The evaluation of the Gaunt factor has been refined 
over recent years since the work of Scheuer (1960) and Os- 
ter (1961) for a range of frequencies, temperatures. Karzas & 
Latter (1961) have given definitive expressions for cjff. Hum- 
mer (1988) has used these expressions to obtain an accurate 
analytical expansion for g// in terms of a two-dimensional 
Chebyshev fit covering a wide range of frequencies and tem- 
peratures. Table|5|lists accurate Gaunt factors (accurate to 
0.7 per cent) derived from Hummer (1988) for a range of 
frequencies and temperatures relevant to CMB studies. 



We now summarize the present situation for the frequen- 
cies (0.4 — 100 GHz) covered by the observations analysed 
here and the observing range of the Low Frequency Instru- 
ment (LFI) on the Planck Surveyor satellite. The free-free 
spectral index is a slow function of frequency and elec- 
tron temperature (Bennett et al. 1992) as shown in Fig. 0| 
Over a wide range in frequencies this can lead to a signif- 
icant discrepancy in the predicted radio emission if this is 
not accounted for. One should therefore still use the accu- 
rate formalism given by Oster (1961) who derives the optical 
depth for free-free emission as 

rc(Oster) = 3.014 x IO'^T^^'^i^ghz" x 
{ln[4.955 x W^iy'^] + 1.51n(Te)} x (£A/)^„,-6pc (6) 

The approximation given by Altenhoff et al. (1960) is often 
used: 

rc(AMWW) = 8.235 x W-'^T-'^-^^u^^^ x (£;M),^-6p, (7) 

For low frequencies (< 1 GHz) the error in this expression 
is of the order a few percent but increases to 5-20 per cent 
for frequencies above 10 GHz. Mezger & Henderson (1967) 
define a factor a{T,i'), the ratio between the two formulae 

as 
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Figure 4. The local free-free spectral index as a function of fre- 
quency and electron temperature; given by equation (3) of Ben- 
nett et al. (1992). The curves from top to bottom are for temper- 
atures Te ={2,4,6,8[heavy], 10,12, 14,16, 18,20) xlO^ K. 



Table 3. The factor a = Tc(Oster)/Tc(AMWW) for relevant fre- 
quencies and temperatures using equation ISl . 



u (GHz) 






(K) 








4000 


6000 


7000 


8000 


10000 


0.4 


0.9972 


0.9934 


0.9912 


0.9889 


0.9844 


1.4 


0.9936 


0.9974 


0.9978 


0.9977 


0.9967 


2.3 


0.9872 


0.9946 


0.9962 


0.9972 


0.9978 


10 


0.9484 


0.9684 


0.9745 


0.9792 


0.9857 


30 


0.8957 


0.9276 


0.9379 


0.9461 


0.9582 


44 


0.8717 


0.9084 


0.9203 


0.9299 


0.9442 


70 


0.8382 


0.8810 


0.8952 


0.9066 


0.9238 


100 


0.8090 


0.8569 


0.8729 


0.8858 


0.9054 



rc(Oster) o.i ^.-0.15 ^ 

{ln[4.995 X 10~Vghz] + l-51n(re)} 



(8) 



Table gives accurate values of a based on equation 
JSJ for a range of temperatures and frequencies. The a fac- 
tor is convenient in that it allows a simple formula to be 
written down using a single spectral index and temperature 
dependence which are modified by the Gaunt factor. 

The brightness temperature Tt is given by 



T41 



TeTc for Tc << 1 



(9) 



For the frequencies relevant to CMB observations and the 
diffuse interstellar medium discussed here, the optically thin 
assumption is valid. At 1 GHz only the brightest Hll regions 
in the Galaxy are optically thick {EM 10® cm~® pc) 
while the Galactic ridge has r ~ 10"^ {EM ^ 10^ cm"® pc 
« 500 R) estimated for Te = 7000 K. 



Finally, the brightness temperature in Kelvin can now 
be written as 

Tb = 8.235 X 10"^are"°-^^!/GHz^(l + 0.08) x {EM) ^^-6p^{10) 

where the factor (1+0.08) is the contribution from the frac- 
tion of He atoms, all of which are assumed to be singly 
ionized; a is given by equation ||HJ. 



3.3 Electron temperatures in Galactic Hll regions 

Optical determinations of electron temperatures in Hll re- 
gions using forbidden line ratios give values in the range 
5000 to 20000 K with a mean value of ?» 8000 K (Reynolds 
1985). These determinations refer to the solar neighbour- 
hood within 1 or 2 kpc of the Sun. 

Electron temperature determinations using extensive 
radio recombination line (RRL) data probe the greater part 
of the Galaxy and show a clear temperature gradient as a 
function of galacto-centric radius. The results of Shaver et 
al. (1983) for 67 Hll regions show a clear gradient in tem- 
perature increasing with galacto-centric radius from 5000 K 
at 4 kpc to 9,000 K at 12 kpc. This is due to the decreasing 
metallicity at larger radii (Panagia 1979) . There is a spread 
of approximately 1000 K at a particular radius. The local 
value at Rq — 8.5 kpc is Te — 7000 K and will be adopted 
as the typical value in the absence of other information. 



3.4 Expected radio emission from Ha maps 

The relationship between radio emission and Hq emission 
can now be calculated using equations Q and 1101 , in units 
of mK/R as 



T 



ff 

_b 



= 8.396 X lO'^a x v. 



-2.lTnO.667i nO-029/T4 

;hz-'4 10 



(1 + 0.08) (11) 



where T4, is the electron temperature in units of 10* K. 
The conversion factors for relevant frequencies are given 
in Table H for Te = 7000 K. Fig. shows (r/^//Hc)i'GHz 
over the frequency range 100 MHz to 100 GHz for Te — 
4000, 7000, 10000 and 15000 K. The increase in (negative) 
slope with frequency is due to the effect of a arising from 
the Gaunt factor dependence. The variation in Tjj-^ at a 
given frequency corresponds to a factor of ~ 2.4 between 
4000 and 15000 K. 

It is of interest to compare our equation IIH with re- 
cent expressions given in the literature. Reynolds & Haffner 
(2002) give tI' /Ihc = TAv^o '^* ^K/i? at Te = 8000 K. 
Our value at 30 GHz and 8000 K would be 17 per cent 
lower at 6.35 jjK/R. This difference may be attributed to 
the 17 per cent over-estimate of the Gaunt factor in equa- 
tion (6) of Valls-Gabaud (1998) which are used by Reynolds 
& Haffner (2002). We have resolved this difference with help 
from David Valls-Gabaud (private communication). The co- 
efficient in his equation (6) should be 3.96 and not 4.4; fur- 
thermore, a typographical error in his equation (11) gave 
the exponent of T4 as 0.317 instead of 0.517. He agrees with 
the calculations of the Gaunt factors in Section r3.2l His new 
calculations also agree well with our equation IIH . Equation 

also from Valls-Gabaud (1998), is not affected. 

The uncertainty in predicting the radio free-free emis- 
sion from the Ha emission as given in equation lllll arises 
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Figure 5. Ratio of radio brightness temperature to Ha intensity 
multiplied by i^q^^^ as a function of frequency. Curves are for 
T =4000,7000[heavy], 10000 and 15000K which is the likely tem- 
perature range corresponding to a factor ~ 2.4 in the ratio. The 
curvature in each curve reflects the frequency dependence of the 
Gaunt factor. 



Table 4. Pree-free brightness temperatures per unit Rayleigh for 
a range of frequencies assuming Tg = 7000 K. 



Frequency 


Free-free brightness temperature 


[GHz] 


per unit Rayleigh 


0.408 


51.2 mK 


1.420 


3.76 mK 


2.326 


1.33 mK 


10 


60.9 AtK 


30 


5.83 /iK 


44 


2.56 tiK 


70 


0.94 imK 


100 


0.43 nK 



mainly from the uncertainty in the electron temperature. 
The ±1000 K uncertainty in the adopted solar neighbour- 
hood value of 7000 K is an upper limit to the scatter of Hll 
region temperatures since it includes some measurement er- 
rors. The corresponding error in equation (IIH is 10 per cent. 
The other possible uncertainty is in whether case A or case 
B applies in the lower density regions at intermediate and 
high Galactic latitudes of interest to CMB studies. For the 
reasons given in Section 13.11 we will adopt the generally 
assumed position that case B is correct and take the asso- 
ciated error in equation I^J to be 1 per cent. However, in 
the unlikely event that case A were to apply, the value of 
Tb/Iua given in equation Ulljl would be increased by 30 per 
cent. Uncertainties in equation Q a-nd in the ionized helium 
to hydrogen ratio in equation llOU are only a few per cent. 
We consider that the overall uncertainty in applying equa- 
tion mil l to be ~ 10 per cent for a large sample of the sky, 
assuming case B applies. 



4 OBSERVATIONAL TESTS OF THE Ha TO 
FREE-FREE RELATION 

We now make an observational test of equation l|ll|l which 
predicts the free- free radio emission for a given Hq intensity. 
In practice, the diffuse free- free emission is mixed with syn- 
chrotron emission so the free-free component must be iden- 
tified and separated. The free-free emission can be identi- 
fied in two ways. Firstly, by its associated radio recombi- 
nation line emission and secondly, by its free-free spectrum 
using multi-frequency observations. Furthermore, a compar- 
ison between the radio and Ha will lead to an estimate of the 
Ha absorption by dust; this can be compared with A{Ha) 
given by equation Q to give a value of the absorbing factor 
of dust fd in front of the Ha. 

4.1 Radio Recombination Lines 

One advantage of using radio recombination lines (RRL's) 
is that the ratio of the line temperature to continuum tem- 
perature can be used to derive an electron temperature 
for the Hll emission which can be substituted in equation 
mil l to give the radio emission expected from the observed 
Ha emission. Two extended regions of Ha emission which 
are sufficiently bright to have been studied in RRL's are 
Barnard's Arc and the Gum Nebula. 

4-1.1 Barnard's Arc 

Barnard's Arc is an ionization-bounded Hll region at 460 pc 
photo-ionized by the Orion I OB star association (Reynolds 
& Ogden 1979). At radio frequencies Barnard's Arc has a 
thermal spectral index (Davies 1963). Gaylard (1984) de- 
scribes RRL observations at 2.272 GHz at 3 positions in 
the Arc. The radio data were converted from antenna to 
brightness temperatures using the factor 1.5 used by Gay- 
lard (1984). The Ha observations from SHASSA and the 
dust map were smoothed to 20 arcmin resolution for this 
study. Table|K|gives the observed brightness temperature Tt, 
Te and Tua, the brightness temperature expected from the 
observed Ha intensity /hq , uncorrected for absorption using 
Te values also given in Table The small (< 10 per cent) 
errors in the radio data for the 3 positions in Barnard's arc, 
indicate that the calculated fd values are accurate to about 
± 0.15 and are therefore consistent with no or little absorp- 
tion [fd ~ 0) in all 3 positions. No absorption would be 
expected if the Ha-emitting gas were nearer than the dust, 
in this nearby (< 500 pc) Gould Belt feature. 



4.1.2 The Gum Nebula 

The Gum Nebula is a diffuse emission region centred on 
{l,b) = (258°, -2°) (Gum 1952). It is a complex at a dis- 
tance of 500 pc containing both free-free and synchrotron 
features which are probably the consequence of supernova 
activity (Duncan et al. 1996; Reynosa & Dubner 1997). The 
RRL data from Woermann et al. (2000) for Ha features 
within the complex taken with 20 arcmin resolution are com- 
pared with the Ha and dust data as in Section [4. 1.11 and are 
given in Table |^ The data tabulated are for the 3 strongest 
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Table 5. Comparison of the observed brightness temperature T;, from RRL data with the predicted brightness temperature Tfjc from 
Ha data using equation llli . The first three positions are in Barnard's Arc (Gaylard 1984) and the next three are in the Gum Nebula 
(Woermann et al. 2000). Using the dust column density D-^ and electron temperature Te, the fraction of dust causing absorption of 
Ha emission is calculated for each observation. 



Coordinates 
(B1950/Galactic) 


Frequency 
[MHz] 


[mK] 


[K] 


[R] 


-^IOOm"! 
[MJy sr-l] 


[mK] 


fd 


5'M9™ +00° 38'" 


2272.661 


295 ± 24 


5200 ± 400 


280 ± 24 


21.2 


330 


-0.1 


5''51™ - 01°00™ 


2272.661 


280 ± 18 


3800 ± 300 


170 ± 15 


15.2 


269 


+0.1 


5''54™ - 02° 45'" 


2272.661 


258 ± 14 


6600 ± 1100 


188 ± 17 


17.2 


254 





G264.00 -12.80 


1715.67 
2422.46 


450 ± 105 
263 ± 53 


8500 ± 2400 
7800 ± 2300 


223 ± 20 
223 ± 20 


6.7 
6.7 


673 
310 


-1.4 
-0.6 


G265.00 -08.00 


1715.67 


570 ± 120 


5700 ± 1400 


232 ± 21 


12.0 


535 


+0.1 


G270.65 +10.40 


1715.67 


398 ± 90 


5800 ± 1400 


123 ± 11 


6.2 


302 


+1.0 



Ha features in the Gum Nebula. There is a large spread in 
fd values from —1.4 to +1.0. The large error bars on the 
radio data {Tb,Te) correspond to an error of approximately 
± 0.6 to ± 1.0 in the derived fd values for the three Gum 
nebula positions. Little can be deduced about the value of 
fd without better radio data; negative values have no phys- 
ical meaning except that the Ha prediction is larger than is 
measured in the radio. 



4.2 Multi-frequency scans across Ha features 

An alternative method of identifying and measuring the ra- 
dio free-free emission from an Ha feature is to use multi- 
frequency scans across it to separate the free-free from the 
underlying synchrotron emission. The brightest features of 
Barnard's Arc are ideal since they lie some 20° from the 
Galactic plane where the background synchrotron emission 
is weak and relatively smooth. We have used the 408 MHz 
map of Haslam et al. (1982) in the destriped form given 
by Davies et al. (1996) and the 2326 MHz map of Jonas 
et al. (1998) to make 3 Galactic longitude scans across the 
Barnard Arc feature at fe = —18°,— 20° and —22° centred 
on I ~ 212°. These are illustrated in Fig. |S| along with the 
corresponding Ha scans from WHAM data; all data are at a 
resolution of 1°. The amplitudes of the best-fitting gaussians 
to the feature centred on each Ha position are given in Table 
|5] Values of fd are estimated as in Section r4.1l for Te = 7000 
K (the average local value) and Te = 5100 K (the weighted 
average of Te from Table |^. The derived fd values are ac- 
curate to about ± 0.1 to ±0.2 and are significantly negative 
when Te = 7000 K is assumed, indicating the radio emis- 
sion is greater than expected from the Ha brightness with 
no foreground absorption. However, when the average mea- 
sured value of Te (5100 K) is used, the derived fd values are 
consistent with no absorption [fd ~ 0). This clearly shows 
the importance of having accurate electron temperatures. In 
this case, the lower electron temperature corresponds to a 
~ 25 per cent reduction in Tua- A further consideration is 
the full-beam to main-beam correction (see Sections 6 & 7) 
which can lead to a factor of ~ 1.5 for objects comparable 
in size to the beam. For Barnard's Arc, which is extended 



(~ 5°) in one direction and ~ 1°.5 in the perpendicular 
direction, the correction is likely to be the square root of 
this factor (~ 1.2). This further correction would lead to a 
slightly positive (~ 0.2) value for fd- 



4.3 Identifying Ha morphology in radio maps 

The morphology of extended bright Ha regions can be iden- 
tified in the low frequency continuum maps even in the pres- 
ence of significant synchrotron emission. By subtracting a 
template from the radio continuum maps proportional to the 
Ha intensity, it is possible to estimate the amplitude of the 
free-free emission. A first-order approach is illustrated for 
the Orion region in Fig.|7| The 408 MHz, 2326 MHz, WHAM 
Ha and the dust maps are first smoothed to 1° resolu- 
tion. Different amplitudes of the free-free emission expected 
from the Ha maps are subtracted from the radio maps until 
the correlation with the Ha maps disappears from the radio 
maps; we are left with the synchrotron features only. Then 
using the absorption expected from the dust map, an 
estimate is made of fd for the field. In this case, assuming 
a nominal electron temperature of Te = 7000 K, both the 
A-Orionis nebula {l,b) ~ (195°,— 11°) and Barnard's Arc 
give fd = —0.2 ± 0.1, suggesting no dust absorption. The 
error bars are estimated by comparing maps with different 
amounts of Ha . As shown in Section 14.21 the effect of a 
lower average temperature for these regions and/or the full- 
beam to main-beam correction, may bring the Ha estimate 
in agreement with the radio data {fd ^ 0). 

A similar study was made of the Gum Nebula where fd 
was found to be —0.2 ± 0.2, and the Ophiuchus Hll region 
I « (0°)-(20°), 6« (+15°) -(+30°) where /d = 0±0.2. The 
increased error in the estimate is due to the complexity of the 
emission across each region. In general, the Ha predictions 
assuming Te — 7000 K are correct when a factor of 0.7 — 1.0 
is included into the Ha estimates (Fig.|7J. 

In summary, the raw Ha predictions in these selected 
regions are slightly too large to fit the current low-frequency 
data by up to ~ 30 per cent. This may be accounted for by 
lower than average electron temperatures (probably the case 
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Figure 6. 408 and 2326 MHz scans across Barnard's Arc at constant b = -18° (left), -20° (middle) and -22° (right). The top panels 
show the 408 MHz scan and the prediction from Ha; both WHAM and SHASSA are shown for each. WHAM is the dotted line and 
SHASSA is the dashed line. The best-fitting gaussian plus a linear baseline is shown as a heavy line; WHAM is the heavier line and 
SHASSA is the lighter line. The middle panels show 2326 MHz data and the same predictions as above for Ha. The bottom panel is the 
scan for D'^ (SFD98). 



Table 6. Comparison of the observed radio brightness temperature Tf, with the prediction from Ha data, Tjio, for 3 scans across 
Barnard's Arc at frequencies of 408 and 2326 MHz. The Ha intensity /jja a-nd dust column density are also listed. The expected 
brightness temperature Thq and the associated dust fraction actually absorbing are calculated for Te = 7000 K (columns 7 &; 8), and 
for Te = 5100 K (columns 9 & 10) (see text). 



Coordinates 


Frequency 


n 




^lOOjim 


Thc [mK] 


fd 


Tea [mK] 


fd 


{l,b) 


[MHz] 


[mK] 


[R] 


[MJy sr-l] 


(Te = 7000 K) 


(Te = 7000 K) 


(Te = 5100 K) 


(Te = 5100 K) 


(212.9°, -18°) 


408 


1980 


67 ±7 


22.5 


3400 ± 300 


-0.6 


2470 ± 220 


-0.2 




2326 


53 


67 ± 7 


22.5 


89 ±8 


-0.5 


64 ±6 


-0.2 


(212.9°, -20°) 


408 


1490 


47 it 5 


30.7 


2400 ± 220 


-0.4 


1730 ± 160 


-0.1 




2326 


53 


47 ±5 


30.7 


62 ±6 


-0.1 


45 ±4 


+0.1 


(212.0°, -22°) 


408 


2990 


80 ±8 


10.8 


4100 ± 370 


-0.7 


2950 ± 270 







2326 


89 


80 ±8 


10.8 


107 ± 10 


-0.4 


85 ±8 


+0.1 
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Figure 7. Morphology of the free-free component in the Orion region, (a) Ha data from WHAM clearly shows Barnard's Arc and 
A-Orionis. (b) 2326 MHz data from Jonas et al. (1998) showing a similar morphology, (c) 100 /jm D'^ map from SFD98 shows some 
correlation but not a one-to-one correlation. Panels (d,e,f) show the 2326 MHz map with varying levels (1.0,0.7 and 0.4 multiplied by the 
conversions given in Table l4l of free- free emission subtracted as predicted by the Ho data. Panels (g,h,i) are the same for the 408 MHz 
map from Haslam et al. (1982). The electron temperature is assumed to be Te = 7000 K. The best subtraction occurs when the factor is 
~ 0.7. Note that no correction has been made for dust in this illustration and that the full-beam brightness temperature correction has 
not been applied (see text). 



in some parts of the Orion region) which will vary the con- 
version factor from Ha to radio continuum by up to ~ 30 
per cent within a temperature range of 5000 — 10000 K. The 
full-beam to main-beam correction may also be a significant 
correction; up to 50 per cent, but ~ 20 — 30 per cent in elon- 
gated structures. The resulting fd values for these regions 
are consistent with little or no absorption {fd ~ 0). This 
might be expected since these regions are at relatively high 
Galactic latitudes (|6| > 10°) and therefore nearby objects 
in which the geometry of the ionization, as seen from the 
Sun's position in the Galactic plane, produces ionization on 
the near side of the gas/dust clouds. 



5 THE BEST-ESTIMATE FREE-FREE 
TEMPLATE 

We are now in a position to move from the observed 
Ha maps of the northern and southern skies to produce 
a full-sky dust-corrected Ha map at an angular resolution 
of 1°. Then, using equation lllll we can derive all-sky free- 
free maps except in areas near the Galactic plane where the 
predicted dust absorption becomes uncertain. 



5.1 Full-sky Ha map corrected for dust absorption 

We obtain a full-sky Ha map by using the northern WHAM 
Ha data at 1° resolution and the higher angular resolu- 
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tion Southern SHASSA Ha data. The Ha map produced 
from WHAM data was made by interpolating the irregularly 
spaced data grid using the Interactive Data Language (IDL) 
routines triangulate and trigrid, which are based on 
Delauney triangulation, as suggested by the WHAM team 
(L.M. HafFner, private communication). The median filtered, 
continuum subtracted composite SHASSA map was used as 
given by the SHASSA team (J.E. Gaustad, private commu- 
nication). Where the maps overlap, the WHAM map is used 
because of its higher sensitivity and better zero-level cer- 
tainty, except for declinations < —15° where SHASSA data 
are preferred. The observed maps are regridded onto an over- 
sampled Cartesian Galactic coordinate grid and smoothed 
to 1° resolution. We then use the map of 100 /im dust 
emission to correct for dust absorption at each point on the 
Ha map according to equation Q. The fraction of the dust 
fd lying in front of the Ha is taken to be 0.33 as shown in 
Section [2.41 and discussed in Section 0] 

Fig. 15.11 shows a colour representation of the dust- 
corrected full-sky Ha map in MoUweide projection. Regions 
where the dust absorption correction is greater than 1.0 mag 
are shown as grey; Ha data become unreliable in these re- 
gions. The Local System (Gould's Belt) is clearly seen be- 
neath the plane at I ~ 180° and above the plane at / 0°, 
reaching to |6| ~ 30° or more. 

It should be remembered that both the WHAM and 
SHASSA surveys are in the early stages of analysis and fur- 
ther revisions are anticipated. Baseline effects are still visible 
in the data as presented in Fig. 15. II For example, the WHAM 
data cannot accurately subtract the geocoronal emission 
near the ecliptic pole {{I, b) ^ (96° , 30°)), while the SHASSA 
data have artifacts due to the varying geocoronal emission 
in each 13° x 13° field which results in a residual power on 
these scales in the region near {{l,b) ~ (320° , —45°)). The 
typical level of these baseline variations are ^ 1 R. 



5.2 A free- free all-sky template 

We can now use Fig. 15.11 to produce an all-sky free-free 
template at any radio frequency using equation llllt . For 
illustration we generate the template for 30 GHz, a fre- 
quency which is widely used for CMB studies such as the 
space projects MAP and Planck and the ground-based ar- 
rays CBI, DASI and VSA. In applying equation IIH we have 
taken Te = 7000 K, the appropriate value for the region of 
Galaxy sampled by the Ha maps at |6| > 10°. At 30 GHz, 
Tt = 5.83 fj.K/R (Table EJl. 

The all-sky template for 30 GHz free-free emission is 
shown in Fig. |U] The amplitude of the free-free signal is 
< 10 nK for l&l > 30° at most Galactic longitudes, although 
in the Local System Tt can be as large as 100 ^K. At higher 
frequencies the free-free emission falls as u~'^'^^ and accord- 
ingly at 70 to 100 GHz, where foreground contamination 
of the CMB is least, the Tb values will be 5-10 times less 
than shown in Fig. |3] At these frequencies free-free and vi- 
brational dust emission are the dominant foregrounds. Only 
the Galactic ridge (|6| < 5° — 10°) and isolated regions such 
as Orion have Tj, > 10 - 20 nK at frequencies 70 - 100 GHz. 

We discuss the contribution of free-free emission to the 
CMB foregrounds in Section f7. 31 



Table 7. Comparison of the r.m.s. fluctuations at 1° resolution 
of 408 MHz data and the free-free template for different Galactic 
latitude cuts. The all-sky value will be a lower limit since the 
corrections for dust are uncertain near the Galactic plane (|fe| < 
5°). 



Galactic 


408 MHz 


Free-free 


Free-free / 


cut 


r.m.s. 


r.m.s. 


synchrotron 




(K) 


(K) 


ratio (%) 


All-sky 


30.3 


2.9 


> 10 


|b| > 10° 


10.3 


0.85 


8 


\b\ > 20° 


7.2 


0.39 


5 


\b\ > 30° 


5.9 


0.23 


4 


\b\ > 40° 


5.1 


0.08 


2 



6 AN IMPROVED LOW-FREQUENCY 
TEMPLATE 

The 408 MHz all-sky map (Haslam et al. 1982) made at 
0°.85 resolution is widely used as a template for the Galac- 
tic synchrotron foreground. Our present study allows us to 
correct this map for free-free emission at 408 MHz to ob- 
tain a clean synchrotron template. The free-free correction 
is estimated from the dust-corrected Ha template of Fig. 15. II 
using Tb — 51.2 mK/R as given by equation IIH and Table 
Hwith Te = 7000 K. 

Over much of the sky the free-free features are weaker 
than the synchrotron features at 408 MHz. However in re- 
gions such as the Local System at intermediate latitudes 
and the Gum Nebula at lower latitudes, the free-free fea- 
tures can be brighter than the synchrotron features. This is 
demonstrated in the Orion region in Fig. |7| where significant 
free-free emission is seen at both 408 MHz and 2326 MHz. 
In such regions the correction is clearly important. 

To estimate the level at which free-free emission is 
present at 408 MHz, the r.m.s. fiuctuation of the 408 MHz 
map smoothed to 1° resolution, are compared with that 
given by the free- free template map. Table [7| gives the re- 
sults for different Galactic latitude cuts. The estimate of 
free-free contamination suggests that at high Galactic lati- 
tude (|6| > 40°), away from bright Ha features, the free-free 
component is a few per cent compared with the synchrotron 
component. At lower latitudes, the fraction increases to ~ 10 
per cent. The all-sky value will be considerably underesti- 
mated since the corrections for dust absorption near the 
Galactic plane (|b| < 5°) are uncertain. These estimates 
agree well with previous best estimates of the ratio of syn- 
chrotron and free-free. For example, the 10,15 and 33 GHz 
Tenerife experiments combined with the 5 GHz Jodrell Bank 
interferometer show that the synchrotron and free-free com- 
ponents at intermediate latitudes are approximately equal at 
~ 10 GHz (Jones et al. 2001). Assuming spectral indices of 
—3.0 and —2.1 for the synchrotron and free- free respectively, 
the ratio at 408 MHz extrapolated from 10 GHz predicts ~ 6 
per cent of 408 MHz emission is in the form of free-free, in 
good agreement with the Ha values given in Table Q 

Higher frequency maps will necessarily have higher frac- 
tions of free-free emission. The widely used maps at 1420 
MHz (Reich & Reich 1988) and 2326 MHz (Jonas 1998) are 
used for estimating the Galactic contribution at higher fre- 




Figure 8. Full-sky dust-corrected Ha map in Mollweide projection. The map has been repixelated into the HEALPix representation (Gorski, Hivon & Wandelt 1999). The colour-scale 
has been histogram-equalized. Longitude / = 0° is in the centre and increasing to the left. Units are log(/?). The composite map uses data from the WHAM survey (North) and the 
SHASSA survey (South) smoothed to a resolution of 1°. Data from WHAM is used in regions where data overlap for declinations > —15°. The data are corrected for dust extinction 
up to a 1 mag of absorption assuming = 0.33 derived in Section 2.4. Absorption above 1 mag is masked off as grey and depicts regions where the true Ha absorption is uncertain. 
Baseline uncertainties are evident at high Galactic latitudes. 
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Figure 9. Free-free brightness temperature template at 30 GHz with 1° resolution. Grey-scale is logarithmic from 5 — 1000 /xK. Regions 
where the template is umeliable are masked white. Contours are given at 5 (dot-dashed), 10,20,40,100,200 and 500 /iK. 



quencies. Using the estimate of 6 per cent at 408 MHz, the 
free-free fluctuations at 1420 MHz will contribute ~ 15 per 
cent compared to the synchrotron component at intermedi- 
ate latitudes, while at 2326 MHz this increases to « 24 per 
cent. By correcting the 408 MHz map for free-free emission, 
a better synchrotron map can be made. This is particularly 
important for cross-correlation analyses which assume that 
the templates are not correlated. 

With Ha data alone it is not possible to construct 
a free-free correction template at lower Galactic latitudes. 
Fig. 15.11 shows regions along the Galactic plane where the 
dust absorption is too large to estimate an accurate dust- 
free Ha intensity. We plan to use multi-frequency maps of 
the Galactic plane to separate the free-free and synchrotron 
components at low latitudes. The narrow free-free latitude 
distribution along the plane has a peak brightness at 408 
MHz which varies from 40 to 80 per cent that of the broader 
synchrotron distribution at I = 10° to 50° (Large, Mathew- 
son & Haslam 1961) . Accordingly, the correction of the 408 
MHz map for free-free emission is important at low latitudes. 

The other factor in constructing a template at 1° res- 
olution from the low frequency radio maps is the effect 
of the main-beam to full-beam ratio (~ 1.5) on angular 
scales of a few degrees; we note that the brightness temper- 
ature scale of the published maps is on the full-beam scale 
5° — 10°). The result is that temperatures on the main- 
beam scale (~ 1°) should be multiplied by a factor ~ 1.5 to 
bring them to the correct brightness temperature scale. This 
procedure has been considered in Sections 14.21 and 14.31 and 
clearly needs to be taken into account for accurate free-free 
and synchrotron predictions. This scale-dependence to true 
brightness temperature has been noted before (e.g. Jonas et 



al. 1998) and is currently under review when applying it to 
existing full-sky maps. 



7 DISCUSSION 

7.1 Precision of the free- free template 

The Ha intensities given by WHAM and SHASSA have 
quoted calibration accuracies of 10 and 9 per cent respec- 
tively. The corresponding zero level uncertainties are ap- 
proximately 0.2 and 0.5 7?. The accuracy is likely to be im- 
proved as these surveys are refined. The largest uncertainty 
in the Ha intensity is in the dust absorption correction as 
discussed in Section |21 At intermediate and higher Galactic 
latitudes where fd x Aua ^ 0.2 mag, the error in the cor- 
rected Ihoi is ^ 5 per cent; this error increases to ~ 30 per 
cent where fd x Ana is 1 mag nearer the Galactic plane. 

On converting the dust-corrected lua to Tt, the free-free 
brightness temperature, a further uncertainty is introduced 
by the spread in Te, the electron temperature, in the local 
region of the Galaxy covered by the Ha maps. For Te = 
7000 ± 1000 K, the Tb uncertainty is ±10 per cent as derived 
from equation lllll . As far as the emission theory used in 
Section|21is concerned, we believe that the Gaunt factor and 
the analytical relations given are accurate to 1-2 per cent. 
It is assumed throughout that the Ha emission is given by 
the case B emission formulae; the observed radio Ti, values 
are consistent with this assumption. 
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Figure 10. Correlation plots between Ha and dust at a resolution of 1° for two regions (a) I = (90° — 120°), b = (20° — 40°) and (b) 
I = (180° — 210°), b = (—20° to —40°). There is an overall positive correlation with considerable scatter. The Orion region (b) is brighter 
and shows a much larger scatter. The Pearson's correlation coefficients are 0.59 and 0.22 respectively. 



7.2 The contribution to other templates 

A radio-independent method of determining a free- free tem- 
plate has a significant benefit in establishing the other 
Galactic foreground templates for CMB studies. We have 
already shown in Section El that an improved synchrotron 
template can be derived from the 408 MHz map by correct- 
ing for the known free-free emission. This method breaks 
down for the small area along the Galactic plane where the 
dust absorption of the Ha becomes too high; other methods 
of determining the free-free emission are available here. 

Spinning dust radiating in the region 10-40 GHz is pro- 
posed as a Galactic foreground (Draine & Lazarian 1998); 
observations appear to confirm this scenario (Kogut et al. 
1996; Leitch et al. 1997; de Oliveira-Costa et al. 1997, 1998, 
1999, 2000, 2002; Finkbeiner et al. 2002). A strong confir- 
mation of the presence of spinning dust radio emission cor- 
related with the SFD98 dust template requires a clean sepa- 
ration of the free-free component beforehand. This is impor- 
tant because the dust and Ha are generally thought of as 
being correlated. An example of a region of the Galaxy where 
the correlation is relatively strong is I — (90° — 120°), 6 = 
(20° - 40°). This is illustrated in Fig.ITOIa) which shows the 
correlation between Ha and (D^) from SFD98 at 1° reso- 
lution. The relation between lua and for the region in 
Orion I = (180° - 210°), fe = (-20° to -40°) is given in Fig. 
llUI (b') where it is seen that the correlation is less strong. The 
distribution of Ha and dust for the Orion region is shown in 
Fig. 111! The situation is clearly complex on the scale of 1° 
as would be expected for a star-forming region and under- 
scores the necessity of removing the free-free emission before 
investigating spinning dust. 



7.3 Free-free contribution to the CMB 
foregrounds 

We will now estimate the free-free contribution to the CMB 
foregrounds at ~ 30 GHz. If not subtracted from the data, 
the foreground signal will add in quadrature to the CMB 
signal in the power spectrum since they will be uncorre- 



lated on the sky. The free-free template allows statistical 
estimates to be made of the free-free component on angular 
scales Z 1° compared to the primordial CMB fluctuations. 
At angular scales of 1° (e « 200), the CMB has r.m.s. fluctu- 
ations of ~ 75 fiK, falling to ~ 30 nK at smaller and larger 
angular scales. Over a large portion of the higher latitude 
sky (|b| > 40°), the average contribution of free-free emis- 
sion is ~ 9 piK at 30 GHz on scales of 1°. This amounts 
to a 1 per cent increase in the power spectrum at ^ = 200. 
However, individual regions may be contaminated at a much 
higher level depending on their position. At higher frequen- 
cies (~ 70 — 100 GHz), the free-free contribution will be 
IS 1 fiK and will be neghgible at high latitudes. 

We can also compare the estimates for specific regions 
of the sky where CMB observations have been made. The 
VSA has observed three regions of the sky measuring angu- 
lar scales of « 0°.3 — 2°.0 {£ « 150 — 900) covering a total area 
of 101 square degrees at a frequency of 34 GHz (Taylor et 
al. 2002). The contamination from free-free emission using 
the free-free template is estimated to give an r.m.s. temper- 
ature of 0.6 — 0.8 at the 1° scale and hence the free-free 
contamination is negligible in these selected regions. 

The North Celestial Pole (NCP) region has been studied 
by the Saskatoon group (WoUack et al. 1997). The Saska- 
toon data show an anomalous component which cannot be 
explained on conventional grounds. Using the free-free tem- 
plate, the r.m.s. variations are ~ 3.4 in a 15°-diameter 
circle centred on the NCP on smoothing the Ha map to 
a resolution of 1°.45 and assuming a conversion factor of 
~ 7 fxK/R. Simonetti et al. (1996) use their Ha data to es- 
timate an upper limit of 4.6 fiK at 27.5 GHz in this region. 
These values of ~ 4/iK can be compared with the « 40 
r.m.s. variations in the Saskatoon data. It therefore seems 
clear, that the anomalous emission is not in the form of free- 
free emission from gas at Te ~ 10000 K. The current view 
is that it may be due to spinning dust emission as proposed 
by Draine & Lazarian (1998). 
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Figure 11. Correlations between Ha and dust in the Orion region. 
3 to 380 R and contoured at 5 (dot-dashed), 10,20,30,40,70,100, 150 
logarithmic grey-scale from 0.1 to 200 MJy sr^-'^ and contoured at 3 ( 
is an overall correlation between Ha and dust but it is not one-to-oi 




(a) Ha data from WHAM shown as a logarithmic grey-scale from 
and 200 R. (b) 100 /xm data from SFD98 depicting the dust in a 
dot-dashed),5,7,10,15,20,25,30,35,40,50,70 and ICQ MJy sr"!. There 



8 CONCLUSIONS 



The recent publication of Ha surveys covering the major- 
ity of the sky has provided a break-through in obtaining 
a free-free Galactic foreground template of importance for 
CMB and ISM studies. Having derived the free-free tem- 
plate it is then possible to determine a better synchrotron 
template and to make real progress in constructing a spin- 
ning dust template. The free-free template has many other 
applications which include a re-analysis of the COBE-DMR 
data (Banday et al., in preparation) and a free- free power 
spectrum analysis (Dickinson et al., in preparation). 



We can look forward to an improved all-sky free-free 
template at higher resolution than the 1° used in the present 
work. A combination of the high resolution (axcmin scale) 
narrow-band filter surveys such as SHASSA and VTSS com- 
bined with the high sensitivity Fabry-Perot survey at a 1° 
scale would be ideal. At lower Galactic latitudes where the 
dust absorption of the Ha emission is significant, multi- 
frequency continuum surveys should be able to separate 
the strong synchrotron and free-free components. Recom- 
bination line surveys such as those being undertaken with 
HIPASS in the South (Barnes et al. 2001) and HIJASS in the 
North (Kilborn 2002) will provide a confirmation of the low 
latitude free-free template and distinguish it from spinning 
dust which will have a very similar latitude distribution. 
The recombination line data will also lead to a kinematic 
2-dimensional picture of the ionized gas distribution in the 
Galaxy using kinematic distances since they contain velocity 
information. 
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